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We investigate adsorption of a gas on the fiat surface of a solid deformable adsorbent taking into account 
thermal fluctuations and analyze in detail the effect of thermal fluctuations on the adsorbent deformation in 
adsorption. The condition for coexistence of two states of a bistable system of adsorbed particles is derived. We 
establish the specific properties of the adsorption-induced surface normal relaxation of an adsorbent caused 
by thermal fluctuations. The mean transition times between two stable states of the bistable system are derived 
in the parabolic approximation and in the general case. 
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1. Introduction 

Investigation of adsorption of gas particles on the surfaces of solids is very important for solving vari¬ 
ous problems of physics and chemistry. Even a monolayer coverage of the adsorbent surface ^ adsorbed 
particles (adparticles) is capable of considerably changing adsorbent properties (see, e.g., 

It is also important to know the amount of adsorbed substance on the surfaces of various adsorbents 
and its distribution over the adsorbent surface for heterogeneous catalysis filllill. 

Various generalizations of the classical Langmuir model of adsorption provide new specific features 
for the amount of adsorbed substance and its kinetics (see, e.g., la llCMl^ l. For example, taking account of 
lateral interactions between adparticles leads to a hysteresis loop of adsorption isotherms and to different 
structural changes of the adsorbent surface I5l-la ll4 - l^ . 

Hysteresis-shaped adsorption isotherms also occur due to the adsorbent deformation in adsorption. 
The effect of adsorption-induced deformation of porous solids is well known [see, e.g., (H) (chapter 4 
and references therein)]. Specific nonmonotonous behavior of the adsorption-induced deformation of 
various porous adsorbents with the gas pressure is observed in a series of experiments; theoretical ex¬ 
planation of this phenomenon and its effect on the adsorption isotherms are given in For a solid 

ideal (nonporous) adsorbent with the flat energetically homogeneous surface, a hysteresis of adsorption 
isotherms of an adsorbed one-component gas caused by the adsorbent deformation in adsorption is es¬ 
tablished in (13 • It should be noted that the hysteresis of adsorption isotherms due to retardation effects 
in the case where the typical adsorption-desoration time is much less than the relaxation time of the 
adsorbent surface was predicted by Zeldovich \ 2 ^ as early as in 1938. 

Memory effects are also essential in the study of the surface diffusion of adparticles over the adsor¬ 
bent surface if the typical time of the moving adparticles is less than the relaxation time of the adsorbent 
surface (see, e.g., (291] and references therein). 

For a bistable system of adparticles on the flat surface of a deformable adsorbent, it is of great interest 
to investigate possible transitions between stable states of a system. One of the ways of correct description 
of these transitions is taking account of fluctuations in a system. This problem is closely connected with 
investigation of the effect of fluctuations on the normal displacement of the adsorbent surface caused by 
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adsorption, especially, taking into account the experimentally established phenomenon of the adsorption- 
induced surface normal relaxation of an adsorbent (see, e.g., the review ( 3 ] and references therein). Only 
a few experimental data for the adsorption-induced normal displacement of the flat surface of an ad¬ 
sorbent are available for some specific values of the surface coverage. According to (23, the theoretical 
dependence of this displacement on the dimensionless gas concentration can be both a continuous func¬ 
tion and a discontinuous function depending on the value of the coupling parameter. The second case 
corresponds to the bistable system under study. 

In the present paper, we generalize the model of adsorption of a gas on the flat surface of a deformable 
solid adsorbent proposed in (^ for the case taking into account thermal fluctuations in the system. In 
section [H we present general relations for the normal displacement of the adsorbent surface in adsorp¬ 
tion and the amount of adsorbed substance on the surface of an adsorbent whose desorption properties 
vary due to the adsorbent deformation and analyze in detail the Important case of the system with two en¬ 
ergetically equivalent states. In section[3l we study the effect of thermal fluctuations on specific features 
of the probability density of the position of the adsorbent surface for monostable and bistable systems. 
The mean transition times of the bistable system between its stable states are investigated in the parabolic 
approximation. In the general case, the corresponding times are obtained in Appendix. 


2. Statement of the problem and general relations 

We consider localized monolayer adsorption of a one-component gas on the flat surface of a solid 
adsorbent within the framework of the model taking into account the adsorbent deformation in adsorp¬ 
tion i^ . Gas particles are adsorbed on adsorption sites with identical adsorption activities located at 
the adsorbent surface. Furthermore, each adsorption site can be bound with one gas particle. The total 
number of adsorption sites N is constant. The Cartesian coordinate system with the Ox-axis directed into 
the adsorbent perpendicularly to its surface is introduced so that the gas environment and the adsorbent 
with clean surface occupy the regions x < 0 and x^O, respectively. 

First we briefly describe the model and the main results in (^ that are necessary in what follows. 

Each adsorption site is simulated by a one-dimensional linear oscillator that oscillates perpendicu¬ 
larly to the adsorbent surface about its equilibrium position (x = 0 in the absence of an adsorbate). Owing 
to the binding of a gas particle with a vacant adsorption site, the spatial distribution of the charge density 
of the adsorption site changes. Furthermore, this change depends on different factors connected with the 
adsorbent and gas particles (see, e.g., (iilii Eil - tiiili^ ). 

As a result, the interaction of the bound adsorption site (the adsorption site occupied by an adparti- 
cle) with the neighboring atoms of the adsorbent on the surface and in the subsurface region changes, 
thus changing the resulting force of the neighboring atoms acting on the bound adsorption site. This can 
be described in terms of an adsorption-induced force F®(r, f), where r is the running coordinate of the 
adsorption site that acts on the bound adsorption site. Due to this force, the equilibrium position of the 
adsorption site (x = 0 for a vacant adsorption site) shifts. Once the adparticle leaves the adsorption site, 
under some conditions, another gas particle can occupy the adsorption site before it relaxes to its nonper- 
turbed equilibrium position x = 0. Thus, a gas particle is adsorbed on the adsorbent surface with changed 
desorption characteristics caused by a local deformation of the adsorbent by the previous adparticle, 
which can be interpreted as adsorption with memory effect. Assume that the adsorption-induced force 
P^{f, t) is normal to the boundary and depends only on the coordinate x: F®(f, t) = F^{x, t) - exF^[x, t), 
where ex is the unit vector along the Ox-axis. 

In the model, the time-step force F^{x,t) acting on the adsorption site during discrete time inter¬ 
vals, where the site is occupied by an adparticle, is replaced by an effective time-continuous adsorption- 
induced force F®®(x, t] - F^{x]6{t] taking into account the presence of an adparticle on the adsorption 
site in the mean. Here, F®(x) = exF^(x) is the adsorption-induced force acting on the adsorption site 
permanently bound to an adparticle and 6(r) = Afb(f)/A( is the surface coverage by the adsorbate and 
A/b(?) is the number of bound adsorption sites at the time t. Expanding F®(x) in the Taylor series in 
the neighborhood of x = 0 and keeping only the first term of the expansion, in terms of the potential, 
F®(x) = -dH®(x)/dx, we have 

V^{x)^-Xx, ( 1 ) 


33603-2 




Adsorption-induced surface normai reiaxation 


where 




dV^ix) 


dx 


x=0 


is the constant adsorption-induced force acting on the bound adsorption site. 

In this mean-field approximation, the kinetics of the surface coverage 8 and the normal displace¬ 
ment X of the plane of adsorption sites, which coincides with the coordinate of a bound adsorption site, 
in localized adsorption, is described by the system of nonlinear differential equations 


dx 

a - \-KX- y6, 

dt 

^ = fcaC(l-0)-fcd(x)0. 
dt 


( 2 ) 


Here, a is the friction coefficient, x is the restoring force constant, C is the concentration of gas particles 
that is kept constant. 



fcd (x) = k- exp 


Edix) 

ksT 


(3) 


are the rate coefficients for adsorption of gas particles and desorption of adparticles, respectively, and 
k- are the preexponential factors, and E^ix) = Ed + xx are the activation energies for adsorption and 
desorption, respectively, Ed is the activation energy for desorption of adparticles from the surface of a 
nondeformable adsorbent (j = 0), T is the absolute temperature, and k^ is the Boltzmann constant. 

The activation energy for desorption Ed{x) depends on the coordinate x due to the adsorbent defor¬ 
mation in adsorption caused by the displacement of the equilibrium positions of bound adsorption sites 
from X = 0. The quantity fcd(x) can be rewritten in the form 


fcdfx) = kd exp 1“ j - (4) 

where the first factor on the right-hand side of l|4) 

fcd = A:_ exp|--^j (S) 

is the classical rate constant for desorption of adparticles in the Langmuir case, which is independent of 
the gas concentration C, and the second factor shows a variation in the desorption characteristic of the 
adsorbent in adsorption of gas particles on its surface. 

In the general case, the activation energy for desonption also depends on the surface coverage 6 due to 
lateral interactions between adparticles (see, e.g., However, even in the absence of lateral 

interactions between adparticles, the model used shows the essential difference of the amount of the 
adsorbed substance on the deformable adsorbent from the classical Langmuir results. 

The first equation of system (2) describes the motion of the bound adsorption site in the overdamped 
approximation ignoring the inertial term of the equation of motion of the oscillator of mass defined both 
by the mass of adsorption site mo and by the mass of adparticle m. This approximation is true if 

« Tr , (6) 

where tm = s/MIx, M - niQ + m, and x^-alx is the relaxation time of an overdamped oscillator. 

The second equation of system 1(2) is the classical Langmuir equation for the kinetics of the surface 
coverage generalized to the case taking into account the adsorbent deformation in adsorption. 

In terms of the dimensionless coordinate of oscillator (or, which is the same, the dimensionless nor¬ 
mal displacement of the plane of adsorption sites) - x/xmax, where x^ax = i/>i; is the maximum station¬ 
ary displacement of the oscillator for the total surface coverage (0 = 1) from its nonperturbed equihbrium 
position X = 0, system m takes the form 

{ (7) 

(10 

— = fcaC(1 - 0) - A:d0 exp (-g^ - 


33603-3 













A.S. Usenko 


where the dimensionless quantity 


rn 

ksT’ 


( 8 ) 


which is called a coupling parameter of adparticles with adsorbent caused by the adsorption-induced 
deformation of the adsorbent or, briefly, a coupling parameter, has the physical meaning of the maximum 
increment of the activation energy for desorption of adparticles (normalized by T) due to the adsorbent 
deformation in adsorption, = l^®(Xmax) = 

Below, we use the model for the case where the variables ^(t) and d{t) are slow and fast, respectively, 
i.e., Tr » tq, where and Tg = TaTd/(Ta -i- Td) are, respectively, the relaxation times of the coordinate 
^(r) and the surface coverage 0(f) in the linear case, Td = l/fcd is the classical Langmuir residence time 
of adparticles (or the typical lifetime of a bound adsorption site) and Ta = l/(fcaC) can be regarded as the 
typical lifetime of a vacant adsorption site. Following the principle of adiabatic elimination (s^] of the 
fast variable 0(f) in 0, we have 

- -r -K’ (9) 

/-Hexp(-gf) 


and the coordinate (f(f) is a solution of the differential equation 


at 

that describes the motion of an overdamped oscillator under the action of the nonlinear force 


( 10 ) 




dUiO 


-< + 


f'-Hexp(-gO 


( 11 ) 


where £ - CK is the dimensionless concentration of gas particles and K- fca/fcd is the adsorption equi¬ 
librium constant in the linear case (x - 0). 

The potential (7(0 can be represented in the form 


( 7(0 = 11 ^( 0 , 


V"(0 = <^-2<--ln 
8 


7' + exp(-g^) 
£+1 


( 12 ) 


Relations 0 and (10) correctly describe the behavior of Of) and 0(f) for times f» Tg. 

Sometimes, instead of equation JTO) . it is more convenient to use the equation of motion of an over¬ 
damped oscillator in terms of the coordinate x. 


dx 


a — = F{x), 
dt 


(13) 


where the force F(x) = - 
lows iiili : 


d(7(x) 

dx 


and its potential (7(x) are expressed in terms of F[^) and (7(0 as fol- 


FU) 

(7(x) 


Vax t^(^) I 


-xx + x 


£ + exp{-bx) ’ 


xx‘- , _, 7' + exp (- b x) 


-r X - Tin - 

2 £+\ 


(14) 

(15) 


where h = x/(A:Br) = g/Xmax- 

In the stationary case, the first equation of system (7) yields 


0 = 0 


(16) 


and the equilibrium position of the oscillator £, describing the stationary displacement of the adsorbent 
surface in adsorption is a solution of the equation 


exp(-g<f). 


(17) 
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According to (i^l, the behavior of system ([7) essentially depends on the values of the control param¬ 
eters £ and g. For g ^ gc = 4, the system is monostable (the single-valued correspondence between the 
concentration £ and the coordinate £, occurs). For g > gc, the system is monostable only for £ $ [£^, 
where £^ and £2 are the bifurcation concentrations defined by the relations 


^fc = (§ffc-l)exp(-g^y. fc=l,2, 

where 

2 ’ ^2-2 

are the bifurcation values of which are two-fold stationary solutions of system Q, the quantity 


(18) 

(19) 


d - 



( 20 ) 


is the width of the interval of instability of the system symmetric about ^ = 1/2, and q - g/4. 

For any £ e {£\, £'^ called the interval of bistabihty, equation (17) has three real solutions < ^2 < fa; 
furthermore, the stationary solutions and ^3 of system (7) are asymptotically stable (i.e., stable nodes) 
while the stationary solution £,2 is unstable (i.e., saddle). 

For g > 4, the adsorption isotherm 6[£) has a hysteresis An example of such an adsorption 
isotherm is shown in figure [T] In this figure and in figure [S] parts of curves of the surface coverage 6{£) 
corresponding to stable and unstable stationary solutions are shown by solid and broken lines, respec¬ 
tively. In view of (16), the curves in figures [l] and also describe fhe displacement £{£) of the adsor¬ 
bent surface from its nonperturbed equilibrium position f = 0 with concentration £, i.e., the adsorption- 
induced surface normal relaxation of a solid adsorbent. 

The surface coverage 6 increases with the concentration £ along the lower stable branch of the 
isotherm ending at the bifurcation concentration £- £2 (here, £^ =; 0.0063 and £2 ~ 0.053); furthermore, 
the increment of 6 depends both on an increase in the gas concentration and on a variation in desorption 
properties of the adsorbent surface caused by the adsorbent deformation. The jump of 6 to the upper 
stable branch of the isotherm at ^ = ^2 is caused exclusively by a change in desorption properties of the 
adsorbent. For convenience, transitions between stable branches of d{£) are shown in figure[T]by light 
vertical straight lines with arrows indicating the direction of transition. Arrows under and above stable 
branches of 6{£) show the direction of variation in £. For £ > £\, the surface coverage 6 varies with £ 
along the upper stable branch. 

The transition of d from the lower stable branch to the upper one at the bifurcation concentration 
is accompanied by a step increase in the activation energy for desorption, which hampers the desorption 
of adparticles from the surface. As a result, as the concentration £ decreases from a value greater than £ 2 , 
the reverse transition of 6 from the upper stable branch to the lower one occurs at the lower bifurcation 
concentration £^ < £2 at which the upper stable branch ends. 



0 0.01 0.03 0.05 I 


Figure 1. Adsorption isotherm for g = 8. 
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This behavior of the surface coverage 0 versus the gas concentration t corresponds to the well-known 
principle of perfect delay 

The essentially different behavior of adsorption isotherms for g ^ 4 and g > 4 depends on the varying 
shape of the function V (0 (also called a potential) in these cases, namely j^: y((f) has a single well for 
g ^ 4 and for g > 4, ^ C [(\, (!\\ and two weUs with local minima at ^ (the first well) and = ^3 (the 
second well) separated by a maximum eXi, - £, 2 , where k- 1,2,3, are the solutions of equation dlT) 
for g > 4, ^ E (^j, ^ 2 )- In the last case, denote Vjt = V{£ic), where k - 1,2,3, and, for a given g, the position 
of the absolute minimum of the potential l^(f) depends on the value of £. In the special case where 
system Q has a two-fold stationary solution for ^ or ^2 for £ - the potential V (f) has a point 
of inflection at ^ A; = 1,2, lying to the right (for fc = 1) or to the left (for A: = 2) of the bottom of the 
single well of l^(f). 

We now dwell in detail on the bistable system with two energetically equivalent states (Vi = V 3 ). For 
any g > 4, this important case occurs for the concentration defined as follows: 

^M = exp(--|j. (21) 

Following (iJliil], the quantity A’m for g > 4 may be called a Maxwell concentration. It is directly verified 
that the potential V{£) with £- £m denoted by 




1 1 cosh2£7^ 

-In-, 

4 2q cosh q 



( 22 ) 


is an even function about £- 1/2. Hence, the function H'^(^) (called a Maxwell potential) at ^ = 1/2 has 
either the minimum (a single well for ^ ^ 1 ) or the maximum (two wells for q>l) 


= (23) 

I 2 I 2q 4 

The double-weU Maxwell potential has equal minimal values denoted by at <5' = ^1 = and 

^ = ^3 = where = (1 + q/q)l2 and 77 is a positive solution of the equation 


X 

— = tanh X. (24) 

q 

This equation directly follows from (17} with £ -£yiOV from the condition that the force F{£) with £ - £m 
denoted by F“(^), 

F“(0 = |(tanh2^f-2f), (2S) 

is equal to zero. According to relations (22) and (^, the wells are separated by a barrier of the height 
Am = - V^. relative to the bottoms of the wells, 

TTi 1 n ’ 


A„=A(i„eosh,-|^). 


2q ■ 


(26) 


In the case where the coupling parameter g is close to the critical gc, i.e., 77 = 1 -t c, 0 < e « 1, we obtain 
77 != \/^ and the very low barrier Am ~ 3 £ 2/8 separating the closely spaced wells. 

For q » 1, we get 77 =; 77(1 -2£u)< which yields ~ 1 - ^m. ~ ^m and 


T/M ~ -/'Z 
‘^min 



In 2 
2q 


(27) 


In this case, the wells are far spaced and the barrier height Am ~ tends to its maximum value equal 
to 1/4 as g —► 00 . 

The curves in figure[2]show the behavior of the Maxwell potential with the coupling parame¬ 
ter g, namely: the growth of the barrier between the wells and the motion of the wells from £- 1 / 2 . 

The curves in figuredepicted for g = 8 and £ e {£^, £ 2 ), where £^ =; 0.0063 and £2 ~ 0.053, clearly 
Illustrate transformations of the double-well potential H(i^) with £. For £ e {£^, £yi], where £m ~ 0.0183, 
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0 0.2 0.4 0.6 0.8 1.0 ^ 



0 0.2 0.4 0.6 0.8 1.0 ^ 


(a) (b) 

Figure 2. (Color online) The Maxwell potential for: (a) g = 4.4 (1), 5 (2), 6 (3); (b) g = 6 (1), 8 (2), 

10 (3). 



Figure 3. (Color online) The potential V((] for g = 8 and the concentration £ = 0.01 (1), £m (2), 0.05 (3). 


the first well is deeper than the second one (curve 1). Thus, the system is stable in the first well and 
metastable in the second one. As£ e i£^, increases, the depths of both wells increase with £ but with 
different increments. For £ - £y^, the double-well potential is symmetric about ^ = 1/2 (curve 2). For 
£ E (/m. ^ 2 ^, the second well is deeper than the first (curve 3). Hence, the system is metastable in the 
first well and stable in the second one. However, according to the principle of perfect delay the 

oscillator, which was initially at rest in the first well, remains in this well with an increase in £ untU 
the well disappears. For the transition of the oscillator from the first weli into the deeper second well 
following the Maxwell principle (^[3, it is necessary to take into account additionai factors that enable 
the oscillator to overcome the barrier between the wells, e.g., fluctuations or the inertia effect (as in (^). 
In the next section of the paper, to investigate transitions of the oscillator from one well into another, we 
take thermal fluctuations in the system into account. 

The bifurcation curve in the plane of control parameters iq, £) in figure|4]divides the first quadrant of 
the piane into two parts. Branches 1 and 2 of the curve correspond to the bifurcation concentrations £^ 
and £ 2 , respectiveiy, defined by reiations (18) and curve M is the Maxwell set of the values of the con¬ 
centration £m for q > 1 defined by (2T)- The domain outside the bifurcation curve is a domain of monos¬ 
tability of the system; for any point of this domain, system (7) has one asymptotically stable stationary 
solution. The open domain enciosed by the bifurcation curve is a domain of bistability of the system; 
for any point of this domain, system (7) has three stationary solutions (two asymptotically stable and 
one unstable). For any point of the bifurcation curve, except for the critical point Pc = (^7c, £c)> where 
qc - I, £c- exp(-2) =; 0.135 is the critical concentration, which is the common point of the branches 
and the cusp of the second kind of the bifurcation curve, system (7) has two stationary solutions (one is 
asymptotically stable and the other one is two-fold). At the cusp, system (7) has one three-fold stationary 
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Figure 4. (Color online) Bifurcation curve (branches 1 and 2 of the curve correspond to the bifurcation 
concentrations and , respectively) and the MaxweU set (curve M). 


solution. Curve M divides the domain of bistability into two subdomains: for any point lying between 
curves 1 and M, the system is stable in the first well and metastable in the second one; for any point lying 
between curves 2 and M, the system is stable in the second well and metastable in the first one. For any 
point of curve M, except for the point Pc, both wells have equal depths and, hence, two stable states of 
the system coexist and curve M is a curve of the coexistence of two stable states of the system. Taking into 
account thermal fluctuations in the system under study, we can conclude that the motion in the plane of 
control parameters [q, £) along any line intersecting curve M is accompanied by the transition of the sys¬ 
tem from one well into another (deeper) well after the intersection of curve M. At the point of intersection 
of the line and curve M, the first-order phase transition occurs j^. 

Given the explicit expression for the Maxwell concentration we can easily plot the actual ad¬ 
sorption isotherms based on the Maxwell principle (and called Maxwell adsorption isotherms) on the 
basis of adsorption isotherms defined by relations (16), (17). To this end, we first note that, for g > 4, 
each stable branch of an adsorption isotherm has a metastable part lying to the right (for the lower sta¬ 
ble branch) or to the left (for the upper stable branch) of the point of the branch with abscissa £ - 
A Maxwell adsorption isotherm consists of the corresponding adsorption isotherm defined by (16), (17) 
without the unstable branch and the metastable parts of the stable branches and the segment AB of the 
vertical straight line £ - £m connecting the stable branches. Hence, a Maxwell adsorption isotherm is an 
adsorption isotherm of a system of adparticles having one stable state for any concentration £ except for 
the sing le concentration £ = £m for which the system has two energy-equivalent states and which, follow¬ 
ing llTll . may be called the phase transition concentration. Note that Maxwell adsorption isotherms thus 
constructed are similar to the well-known adsorption isotherms taking into account attractive lateral in¬ 
teractions between adparticles on a nondeformable adsorbent (see, e.g., islR Iliill^ ). Furthermore, the 
Maxwell concentration fHS agrees with the corresponding expression for the phase transition pressure 
for the classical Fowler-Guggenheim adsorption isotherms in the case of attractive lateral interactions of 
adparticles on a nondeformable adsorbent if the maximum increment of the activation energy for 
desorption of adparticles \V^\ caused by the adsorbent deformation in (8) is replaced by the modulus 
of the energy of attractive lateral interaction of two neighboring adparticles for 6 = 1 multiplied by the 
coordination number of the lattice of adsorption sites. 

The adsorption isotherm and the MaxweU adsorption isotherm consisting of stable branches of the 
adsorption isotherm without metastable parts and the segment AB of the vertical straight line £ - £m 
connecting the branches in figure [S] clearly illustrate the essential difference (the presence and absence 
of a hysteresis loop, respectively) in the behavior of the amount of adsorbed substance with £. The vertical 
part (the segment AB) of the Maxwell adsorption isotherm indicates the coexistence of two stable states of 
the system under study and corresponds to the first-order phase transition in the system. It is easy to see 
that the areas of the domains enclosed by the adsorption isotherm and the segment AB to the left and to 
the right of the segment AB are different. Nevertheless, the equality of these areas is proved. Furthermore, 
for any g > 4 (and, hence, the well-known MaxweU rule of equal areas is true for the adsorption isotherms 
on a deformable adsorbent) if Inf is laid off along the abscissa axis instead of £ (see also the MaxweU rule 
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Figure 5. Adsorption isotherm (thin curve) and the Maxwell adsorption isotherm (thick curve) for g = 8. 

of equal areas for the adsorption isotherms taking into account attractive lateral interactions between 
adparticles on a nondeformable adsorbent (l^). It is also worth noting that the ordinate of the point of 
intersection of the unstable branch of the adsorption isotherm and the vertical straight hne ( - 
equal to 1/2 for any value of g due to the evenness of the Maxwell potential ^"^(0 about ^ = 1/2. 

According to the remark before figure[T] the thick curve in figurej^also describes the actual (Maxwell) 
behavior of the adsorption-induced surface normal relaxation of a solid adsorbent with concentration £. 


3. Probability density 


To investigate transitions of the bistable system between its stable states with variation in the control 
parameters £ and g, we take into account thermal fluctuations in the system introducing a Langevin 
force .^(r) in the right-hand side of the deterministic equation of motion of oscillator IsTH^ which 

yields the stochastic differential equation 


a— - F{x) + ^[t). 

at 


(28) 


The random force ^{t) has the properties of a white noise: 


(.^(f)> = 0, {^[t)^[t'))^2kjiTa8[t-1'), 


(29) 


where the angular brackets (...) denote the averaging over an ensemble of realizations of the random 
force the quantity 2k^Ta in the correlation function in is the intensity of the Langevin force, 
and 5{x) is the Dirac 5-function. 

Following ii^ . we denote random variables by capital letters and their values by small letters (for ex¬ 
ample, X(r) is a realization of the dynamical variable x at the time t). Equation |28) can be reduced to the 
Fokker-Planck equation for the probability density p[x, t) - (5 [x-X{t])} of the coordinate of oscillator 
which also describes the adsorption-induced surface normal relaxation of a solid adsorbent with 
regard for thermal fluctuations. 


a 


dp{x, t) 
dt 


d 

dx 


dU[x) 

dx 


p[x, t) + k^T 


dp[x, t) 
dx 


(30) 


Given the quantity p{x, t) as a solution of the Fokker-Planck equation <30t . the probability den¬ 
sity p{£,, t) of the random variable E is expressed in terms of p(x, t) as follows: 


p(f, t) = IXmaxI p{x, t) I (31) 

By virtue of (9), the joint probability density p(^,6; t) of the random variables and H and 0 has the 
form 

p(^,0;t)^p{^,t)S[0-f{O), 0e[O,l], (32) 
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where the 5-function on the right-hand side of is the conditional probability density 

/(O) with the sharp value for 6 = f{^] and f{^] Is the deterministic function equal to the right-hand side 

of relation 

We first consider the stationary case. Under the natural boundary conditions, the stationary proba- 
blLity density p[x) has the Boltzmann distribution 


p{x)-jY exp 


U{x) 


knT 


(33) 


where J{ is the normalization constant defined as follows; 


OC 


jV - j dxexp 


U{x) 


k^T 


(34) 


In view of f33) . (12) . and fTS) . we get 


p((f) = ^exp|--V'(Oj, 


where 


^ — I I ^ — 


g 1 + t 




(33) 


(36) 


271 1 + A’ 

According to ([33), the functions p{0 and V(^] have extrema at the same points, moreover, if U(0 has 
a minimum (maximum) at some point, then p[^) has a maximum (minimum) at this point Hence, 

the stationary probability density p{^) is single-modal for g ^ 4 and for g > 4, ^ C [(\, and blmodal 
for g > 4, ^ e (fj, (\). 

By using the explicit expression for V{^), we obtain 


piO- 


exp 


ge \ l + /exp(gO 
2 1 1-tA 


\ + i exp [bx] 

pix) =-^— - Poix), 

1 + A 


(37) 


where 


Poix) = 


\/27T CTo 


exp 


,-2 


2 a: 


ol 


(38) 


is the Gaussian distribution of the probabUity density for a linear oscillator with zero mean ((X> = 0) and 
the variance a^- k^Tlx. 

We first consider the single-modal stationary probability density p{x). In this case, the random vari¬ 
able X has the nonzero mean 

(X> = -^x„,ax, (39) 

l-t A 

which yields sign(X> = sign;t; and, hence, the maximum of the probability density pix) shifts in the 
direction of the action of the adsorption-induced force; the variance a^ = iiX - (X»^>, 


a 


1 + g 


(1 + A )2 


(Tn 


(40) 


which is greater than cTg for any values of the concentration and the coupling parameter and, for a fixed 
value of g, reaches its maximum value equal to (1 -i- q)a^ for A = 1; the asymmetry ratio S = {iX - 


(X»3>/ct2, 


^ 3/2 A(l-A2)(l-rA)2 . 

, —TTT X, 


(41) 


[(l + A)2 + gA]3 

and, hence, sign S = sign ((^m “ ^)x)> which implies the change in the sign of the asymmetry ratio in 
crossing the concentration £m', and the excess (flatness) E - ({X - (X»^>/cr^ - 3, 


E^g' 


2 A(H-A2-4A) 


[(l + A)2 + gA] 


2 • 


(42) 
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According to (42), the probability density p[x) is flat-topped (E < 0) for ( e {(-,(+), where C± - 
(2 + 1 / 3 )^M. or peaked (E > 0) for £ t [(-,£+] relative to the Gaussian distribution with mean and 
variance d^ . For ( - (- and ^ , f? = 0 as for this Gaussian distribution. 

In the Maxwell case {£ - ^m). the probability density p{x) in denoted by pu^x) is simplified to the 

form 


Pm{x) = 


/'I /4 


\/271 (Tq 


exp 


2 ^0, 


cosh 






(43) 


which is an even function about x - Xmax/2. 

In the single-modal case (g ^ 4), we have (X) - Xmax/2, + S= 0, audit = -2r7^/(l-i-^)^ < 0. 

Hence, puM is a flat-topped distribution symmetric about its maximum value equal to £^^/[\/2n(jQ] 
at X = Xmax/2. 

We now investigate the transition of the bistable system from one stable state to another due to ther¬ 
mal fluctuations. This transition occurs for the double-well potential V{^] studied above: the left-hand 
and right-hand weUs with minimal values Vl = Vi and (4 = ^ 3 . respectively, at ^ and f = fa = f r 

are separated by the barrier with maximum value 14 = V 2 at f = f 2 = fs • By using the Gardiner repre¬ 
sentation of the Kramers method a , the following system of equations is derived from the Fokker- 
Planck equation J30) : 


dPrit) 

— -Tz — - -khPhW + k^PB.{t), 
at 

dPRit) 

-^—^-k^PM + k^PM. 


(44) 


This system describes the evolution of the probabilities of the presence of an oscillator to the left (in the 
left-hand well), P^it), and to the right (in the right-hand well), P^it], of the point fe at time t under the 
assumption that the relaxation times of the oscillator in the wells are much less than the mean transition 
times between the wells 




00 

PL(f)= / dfp(f,r), PR(r) = Jdfp(f,r), PM + PM^l. 


Here, 




1 


i = L,R, 


(43) 


(46) 


gpPitr ’ 

is the coefficient of the escape rate of an oscillator from the ith well (f = L, R) derived under the assump¬ 
tion of a vanishing probability of the presence of the oscillator in the well outside a small neighborhood 
of f i, f = L, R, and Pi is the probability of the presence of the oscillator in the ith weU (i -L, R) in the 
stationary case, i.e.. 




Pl= dfp(f), Pr= dfp(f) 




and 


- 

J P(f) 




The solution of system (|44) has the form 


Pi(f) = Pyexpl-—l+Pi 


l-exp|- — 


i = L, R, 

where P9 is the initial value of the probability Pi R, for f = 0, A;+ = A;l - 1 - fcR, and 


(47) 


(48) 


(49) 


r+ = l/fc+ = gpPLPRJr 


(SO) 


33603-11 














A.S. Usenko 


is the relaxation time of the quantity i - L, R. 

According to (44), the mean transition time from one stable state at ^ f, to another at , where 
i,j - L, R, i + j, is defined as Ti^j - 1/ki and, with regard for (46), has the form 


Ti^j^g^PiTi, i ,7 = L,R, i'A;, 


(SI) 


which yields the well-known relationship 0 between the probabilities P, and the mean transition 
times 


Pl 

^’r 7r—l 


(S2) 


In the parabolic approximation of the potential U (^) (12) in the neighborhoods of its extrema, the 

mean transition times Ti^j are estimated as follows: 


where 


2ti 


^expl^Afi,/], i,j = L,R, / A 7 , 


A i,j = Vi - Vj = f- 2) + - In. i ,7 = L, R, B, 

g 7 

L(f,g) = l + ga^-l) 


is an even function about ^ = 1/2 such that, for g> 4,1 > L[^i,g) >0,i- L,R, and L(^B,g) < 0. 
In this approximation, the quantity 


(S3) 


(54) 

(55) 


Tr 

Uii,g)’ 


/ = L,R, 


(S6) 


has the sense of the mean relaxation time of an overdamped oscillator in the left-hand (/ = L) or in 
the right-hand (i - R) parabolic potential weU centered at if = and derived from (12). Hence, in this 
approximation, the restoring force constant x is simply replaced by ;><,■ = xL{f,i,g), i - L,R, that already 
depends on d and g, which yields t/ > Tr, / = L, R. 

If the concentration ( d^) approaches the end point of the interval of bistabihty (d^ or d^), then 

i('fRig) or i(^Lig). respectively, tends to zero. According to (^, this is accompanied by an essential 
increase in the relaxation time tr (in the first case) or tr (in the second case) due to the flattening of 
the corresponding well. This indicates that the parabolic approximation is false in small neighborhoods 
of the bifurcation concentrations d^ and d^ because the very shallow right-hand (in the first case) and 
left-hand (in the second case) wells are almost flat-bottomed. 

Formally introducing the quantity 


Tb 


Tr 

|L(fB,g)l’ 


(37) 


which can be regarded as the mean relaxation time of an overdamped oscillator in the parabolic po¬ 
tential well centered at f = ^b and derived from (12) with |L(fB>g)l instead of L(fB,g). we can express 
the mean transition times Ti^j (5^ in terms of the relaxation times of an overdamped oscillator in the 
corresponding parabolic potential wells as follows: 


TiX 271 s/Ti Tb exp|^ Ab,, /,7 = L,R, / A 7 , (38) 

which agrees with the classical Kramers formula in the overdamped case 0. 

The assumption of two time scales (the short time scale for the relaxation of the oscillator in the 
weU where it was at the initial time and at the long-time scale for the transition of the oscillator from 
one stable state to another across the unstable state) used in the derivation of system (44) imposes the 
foUowing condition on the Arrhenius factor: 


exp 


(fA..,) 


» 


271 


(39) 
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Substituting in and JS 8 ), we express the ratio of the mean transition times between the wells 
in terms of the coordinates of their minima and as follows: 

expff Ar,l)= r/^|^exp[f (^R-^L)(fR + <L- 2 )]. (60) 

For the Maxwell potential (22), we get (fs = 1/2, the quantity L(^B,g) = 1 “ ^ is independent of the 
concentration and, hence, tb = T^Kq- 1 ), 



i(fL,g)=i(^R,g)=iM(g)= l-t? l-\ . Tl = Tr=Tm= (61) 

I q‘^1 LM(g) 

where, as above, q is a positive solution of equation (TM . equal mean transition times between the wells, 
Tl^r = Tr^l = 7m and, with regard for (54) . and (63, 

j q^\ 

Tm ~ 2n x/tm Tb coshq exp-. (62) 

V 2ql 

It is worth noting that expression (63 cannot be used if q is close to 1 because, in this case, the 
assumptions of derivation of Fm are violated. 

For large values of the coupling parameter, q » 1, the quantity 77 =: q (1 - 2 which yields Tm ~ Tr, 

Tb ~ Tr/q, 

FM-Tr^expf^]. (63) 

In terms of the potential U{x), we have U{x- Xmax/2)/fcB T =; ql2 for q » 1, which leads to the exponen¬ 
tial growth of the mean transition time (63 with an increase in the coupling parameter g. 

Relations for the mean transition times derived in the general case are given in Appendix. 


4. Conclusions 

In the present paper, we have investigated the problem of adsorption of a gas on the flat surface of 
a solid deformable nonporous adsorbent with regard for thermal fluctuations and analyzed the effect of 
thermal fluctuations on the normal displacement of the adsorbent surface and, hence, on the amount of 
adsorbed substance. We have derived exphcit expressions for the mean transitions times between stable 
states of the bimodal system and investigated the dependence of these times on the values of the coupling 
parameter and the gas concentration. 

According to the estabhshed results, the behavior of the system under study is most interesting for 
the case where the system is bistable. For the bistability of a system, the coupling parameter must exceed 
the threshold value. Thus, first of all, the value of the coupling parameter for the investigated adsorbent- 
adsorbate system should be calculated. However, this requires an additional information because the 
coupling parameter is expressed in terms of the phenomenological constant adsorption-induced force j. 
Nevertheless, this unknown force (and, hence, the coupling parameter) can be expressed in terms of the 
maximum change in the first interplanar spacing x^ax in the case of the total monolayer coverage of the 
adsorbent surface by adparticles: 

X = g^XX^^/kjiT. (64) 

Having the experimental value of Xmax for such an adsorbent-adsorbate system measured with proper 
accuracy, the value of the coupling parameter is easily calculated by using the second relation in (M) . 
Given a base of experimental data (lacking at present) for Xmax for various adsorbent-adsorbate sys¬ 
tems, one can select systems for which the established effects of bistability of the system caused by the 
adsorption-induced deformation of the adsorbent are possible. 

Since the coupling parameter is proportion to x^^^ J64l . experiments with sohd nonporous adsorbents 
with flat surface should be performed for materials with a considerable value of the adsorption-induced 
normal displacement of the surface. 


33603-13 







A.S. Usenko 


Note that the results established in the present paper have been obtained for the model imposing the 
restrictions on the values of characteristic times (the average time between collisions of gas particles with 
an adsorption site, the average residence time of an adparticle on the surface, and the relaxation time of 
a bound adsorption site) and the friction coefficient. 

Since the proposed model of adsorption on a deformable solid adsorbent does not take into account 
lateral interactions between adparticles, it is of interest to generalize this model to the case taking into 
account the joint action of both factors (adsorption-induced deformation of the adsorbent and lateral 
interactions between adparticles). 
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Appendix 


The mean transition times Tj^j derived above in the parabolic approximation can also be deter¬ 
mined for the general case. By using expression (12) for V ((f) and relations (^-(37), ll46)-(48). and l l^ . 
we obtain the following relations for stationary probabilities: 

2 (i + A) [erfc(v^^B) + Aerfc(v^(<B-l))], Pl^I-Pr, (A.l) 

and the mean transition times 


Ti^j^T,Pi2[l + A.) 


=0 


exp 




exp 


Here, 


{2^7(fL-fR)[2(u+l)-fL-^R]}i"(\/^(?i+l-fR))}. Gi = L,R, 

oo 

erfc(z) = dt exp (- 1^) 

\7l J 


(A.2) 


is the complementary error function the special function 


Z. 

A'(z) = exp (-z^) Jdtexp[t^), 


which is expressed in terms of the error function of imaginary argument, is bounded for all real z, has 
the maximum value ;= 0.541 for z 0.924, the known power expansion for |z| < oo, and the asymptotic 
formula F[z] ~ l/( 2 z) as z —► oo. 

In the Maxwell case (f = ^m). relations dA.D yield equal values of the probabilities of the presence of 
the oscillator in both wells: Pl = ^’r = 1/2, which is quite natural for the symmetric Maxwell double-weU 
potential (22). Relations dA.2) are reduced to the Maxwell mean transition time 

Tm = Tr4 v^expf-^] ^ (A.3) 

^ ^ ^ n=0 


where 

Sn (^?) = P(/3nJ - exp [-7 (2 - 7 )] P(^„ - 7 ), 

77 is a positive solution of equation (24). 


Pn ■ 


■(2n+l), 



(A.4) 
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For q » 1, relation <A.3t is sintplified to the form 

(q\ 

rM=:Tr--^exp|^J. (A.5) 

Comparing lA.St and f63) derived for the general case and in the parabolic approximation, respec¬ 
tively, we see that, for large values of q, the mean transition times differ only by the factor \/nl2. 
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iHAyKoeaHa aACopGi^ieio HopMa,nbHa pe/iaKcai^iji noeepxHi 
TBepAoro aACop6eHTy 

O.C. VceHKO 

iHCTHTyr TeopeTMHHOi c|)i3MKH iM. M.M. Boro.nKD6oBa HAH VKpaiHi/i, 

Byn. MerponotHHa 14 6, 03680 Kmib, YKpaiHa 

flocniflxeHO aflcopbitiK) rasy Ha nnocKiii noBepxHi TBepfloro aflcop6eHTy, akhh fletfiopiviyeTbcti npi/i aflcopOpM, 
BpaxoBytOHH rennoBi (|)nyKTyaL(M, i fleranbHO npoaHanisoBaHO Bn./ii/iB Ten.noBHx (|)nyKTyaL(iH Ha flec|)opMaLtiKD afl- 
copOenty. OrpHMano yMOBy cniBicnyBaHHB flBOx craniB 6iCTa6i.nbH0T CHCreMi/i aflCop6oBaHHx nacTi/iHOK. Bcra- 
HOB.neHO oco6nnBOCTi inflyKOBanoK) aflCopOpieio HopManbHOi pe.naKcaLtiT noBepxni aflCopbeHiy, tixi o6yMOB.ne- 
Hi Ten.noBHMH c|).nyKTyaLtitiMH. SnaiifleHO cepeflHi nacH nepexofliB Mix flBOMa ctIhkhmh CTanaMi/i 6iCTa6i.nbH0T 
CHCreMH B napa6o.niHHOMy Ha6.nnxeHi i b 3ara.nbHOMy BwnaflKy. 

KnioMOBi cnoBa: aficopSi^ift, iaorepMa, flecpopMapio, (pjiyKTyapil 6icTa6inbHicTb, ricrepeai^c 


33603-16 





